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Abstract 



O 

O. 

, Context. Timing analysis of Accretion-powered Millisecond Pulsars (AMPs) is a powerful tool to probe the physics of compact 

^v^j ■ objects. The recently discovered IGR J175 1 1-3057 is the 12 discovered out of the 13 AMPs known. The Rossi XTE satellite provided 

an extensive coverage of the 25 days-long observation of the source outburst. 

Aims. Our goal is to investigate the complex interaction between the neutron star magnetic field and the accretion disk, determining 
r_jj , the angular momentum exchange between them. The presence of a millisecond coherent flux modulation allows us to investigate such 

interaction from the study of pulse arrival times. In order to separate the neutron star proper spin frequency variations from other 
HH ■ effects, a precise set of orbital ephemeris is mandatory. 

Methods. Using timing techniques, we analysed the pulse phase delays fitting differential corrections to the orbital parameters. To 
remove the effects of pulse phase fluctuations we applied the timing technique already successfully applied to the case of an another 
AMP, XTE J 1807-294. 

Results. We report a precise set of orbital ephemeris. We demonstrate that the companion star is a main sequence star. We find 
pulse phase delays fluctuations on the first harmonic with a characteristic amplitude of about 0.05, similar to what also observed 
in the case of the AMP XTE J1814— 338. For the second time an AMP shows a third harmonic detected during the entire outburst. 
The first harmonic phase delays show a puzzling behaviour, while the second harmonic phase delays show a clear spin-up. Also the 
third harmonic shows a spin-up, although not highly significant (3cr c.L). The presence of a fourth harmonic is also reported. In the 
hypothesis that the second harmonic is a good tracer of the spin frequency of the neutron star, we find a mean spin frequency derivative 
for this source of 1.65(18) x 10~ 13 Hz s" 1 . 

Conclusions. In order to interpret the pulse phase delays of the four harmonics, we applied the disk threading model, but we obtained 
different and not compatible M estimates for each harmonic. In particular, the phase delays of the first harmonic are heavily affected 
by phase noise, and consequently, from these data, it is not possible to derive a reliable estimate of M. The second harmonic gives a 
M consistent with the flux assuming that the source is at a distance of 6.3 kpc. The third harmonic gives a lower M value, with respect 
to the first and second harmonic, and this would reduce the distance estimate to 3.6 kpc. 
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^ ; 1. Introduction an uncertainty of 0.6". IGR J175 11-3057 was observed by Swift 

• i— i , . (Bozzo et al. 2009), producing a description of the X-ray spec- 

, Accretion-powered millisecond pulsars (hereafter AMPs) are ^ TolTef , gt aJ (2009) ^ a Me near . infrared coun . 

"£J ■ transient low mass X-ray binaries, which show a coherent mod- terpart A deta0ed spectral analyfjis and a set of orbkal parame . 

Ctf ; ulation of their X-ray fluxes with periods of the order of few terfj were iyen b p ift0 et al (201Q) anal m a XoQ XMM _ 

milliseconds. In the recycling scenario AMPs are seen as the Newtm observation . Ri io et al . (2009)? ana i ysing a RXTE ob _ 

progenitors of the millisecond radio pulsars (see e.g. van den servation ^ refined the orbital parameters . Very recently Miller- 

Heuvel 1984), the accretion process being responsible for the Jonef , gt aJ (2Q09) ^ an ^ Qn ^ radio emission 

spinning up of the neutron star (hereafter NS) to milliseconds Surprisingly? another bmsient x . ray source (XT E J1751-305) 

penods - went into outburst very near the position of IGR J1751 1-3057 

The AMP IGR J1751 1-3057 was discovered by INTEGRAL Qn ? Qctoberi 2009 (C henevez et al. 2009) and its pulsations 

on 12 September, 2009 during a galactic bulge monitor- was dete cted by /JXre (Markwardt et al. 2009a) while observing 

ing (Baldovin 2009). Although very close to the previously i GR ji 75 i i_ 305 7. IGR J1751 1-3057 has faded under detection 

known AMP XTE J1751-305, the source position measured by threshold on g October 2009 

INTEGRAL suggested it was a newly discovered X-ray source. ^ ^ WQrk wg resent a detailed timi anal is of me 

The observation of a coherent modulation of the X-ray flux in RXTE ToO observation of the source IGR J1751 1-3057. 
the data from a ToO observation performed by the Rossi X- 
ray Timing Explorer (hereafter RXTE) with a period of about 4 

ms (Markwardt et al. 2009b) permitted to classify IGR J1751 1- 2 Observation and Data Analysis 

3057 as an AMP and confirmed it as a new transient X-ray 

source. Altamirano et al. (2010) reported the presence of burst In this work we analyse RXTE observation of IGR J175 1 1-3057. 

oscillations at the NS frequency. An analysis of a Chandra obser- In particular, we use data from the PCA (proportional counter 

vation by Nowak et al. (2009) gave the best source position with array) instrument on board of the RXTE satellite (Obsld P94041 
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Table 1. Orbital and Spin Parameters for IGR J175 11-3057. 



5 10 15 20 25 30 35 40 
Time (MJD - 55088.0) 

Figure 1. PCU 2 count rate (2-25 keV), subtracted of the background, 
is reported as a function of time during the outburst. The superim- 
posed model represents the best fit using a piecewise linear function. 
The abrupt flux raise at the end of the outburst is due to the onset of 
an outburst from the AMP XTE J1751-305. The dashed vertical lines 
are in correspondence of the type-I bursts present in the observation 
while the continuous vertical lines are in correspondence of the slope 
change accordingly to the model used to describe the count rate. The 
filled circles are relative to the Obsld 94041, while the triangles are the 
XTE J175 1-305 observation (Obsld 94042) 



and P94042). We used data collected in event packing mode, 
with time and energy resolution of 122fis and 64 energy chan- 
nels respectively. We selected data in the energy range 2-25 keV 
in order to maximise the signal to noise ratio, since above ~ 20 
keV the background dominates. The X-ray flux follows a piece- 
wise linear decay as showed in Figure 1, with a peak flux of 
70 counts s _1 PCU 1 . The re-brightening visible in Figure 1, 
24 days after the start of the observation, is due to the AMP 
XTE J175 1-305 in the field of view of RXTE going into out- 
burst. The analysed data cover the time span from 12 September 
2009 (MJD 55086.8) to 22 October 2009 (MJD 55126.8). We 
corrected the photon arrival times for the motion of the Earth- 
spacecraft system with respect to the Solar System barycentre 
and reported them to barycentric dynamical times at the Solar 
System barycentre using the faxbary tool (DE-405 solar system 
ephemeris). We used the Chandra source position reported by 
Nowak et al. (2009), and reported in Table 1 . The uncertainty on 
the source position quoted by Nowak et al. (2009) is 0.6", Icr 
confidence level. 



2.1. Derivation of the orbital ephemeris 

To obtain a first estimate of the mean spin frequency we con- 
structed a Fourier power density spectra of the first data file of 
the ObsID 94041-01-01-00 with a bin size of 2" 11 seconds and 
on data segments of 64 seconds and averaging 53 power spectra. 
We found a strong signal at ~ 244.81 Hz, in good agreement 
with the value reported by Markwardt et al. (2009b). We divided 
the observation in time intervals of about 400 seconds each and 
performed an epoch folding search on each data interval around 



Parameter 


Value 


RA (J2000), (Nowak et al. 2009) 


17 h 51 m 08 s :66 


Dec (J2000), (Nowak et al. 2009) 


-30° 57'41'.'0 


Orbital period, P or t (s) 


12487.5121(4) 


Projected semi-major axis, a x sin i (lt-ms) 


275.1952(18) 


Ascending node passage, T* (MJD) 


55088.0320279(4) 


Eccentricity, e 


< 3 x 10~ 5 


Mass function ', f x (M Q ) 


1.070854(21) x 10~ 3 


Reference epoch, T (MJD) 


55088.0 


Mean spin frequency, v (Hz) 


244.83395156(7) 


Constant v Model best fit parameters 




^0r 2 /d.o.f.) 


1.74(238.5/137) 


Spin frequency at T , vo (Hz) 


244.83395145(9) 


Spin frequency derivative, v (Hz s~') 


1.45(16) x 10~ 13 


Physical Model best fit parameters 




^Ct 2 /d.o.f.) 


1.70(232.8/137) 


Spin frequency, v (Hz) 


244.83395145(9) 


Accretion rate at T ( M Q year 1 ) 


0.92(10) x 10" 9 



Errors are intended to be at lcr c.L, upper limits are given at 95% 
confidence level. Times are referred to the barycentre of the Solar 
System (TDB). Best fit spin parameters are derived in both hypothesis 
of a constant spin-up and flux dependent spin-up, and the uncertainties 
on the given values of v, v and M max include systematics due to the 
uncertainties in the source position (see text). Here we report only the 
second harmonic best fit spin parameters. 

(1) This value was obtained using the latest available measure of G, c 
(http://physics.nist.gov/cuu/Constants/) andM Q (http:// 
nssdc . gsf c . nasa . gov/planetary/f actsheet/sunf act . html). 



the averaged spin period with a period resolution of 4 x 10" 9 s. 
For each time interval we obtained an estimate of the best spin 
period. We excluded all the intervals for which the maximum in 
the^f 2 curve was not significant (3cr c.L), according to the crite- 
rion stated by Leahy et al. (1983). A sinusoidal Doppler modu- 
lation of the spin period due to the source motion in the binary 
system was evident. We fitted the Doppler frequency shifts with 
the formula: 

v(t) = v + v(t- T ) C os l(t), (1) 

'orb 

where vo is the spin frequency at the time To, v is the spin fre- 
quency derivative, A is the orbit projected semi-major axis over 
the speed of light and /(f) = 27r(f-r*)/P or b, where is the time 
of passage through the ascending node and P orb is the orbital pe- 
riod. With a reduced^ 2 (hereafter^ 2 and defined as^ 2 / d.o.f.) of 
0.59(532.2/898) we obtained a first set of orbital parameters and 
a much better estimate of the barycentric spin frequency. 

Using this preliminary orbital solution we analyse the pulse 
phase delays to get a more precise estimate of the orbital and 
spin parameters. We epoch folded data on time intervals of about 
1500 s using 32 phase channels. An example of the folded pulse 
profile is reported in Figure 2. An harmonic decomposition of 
each pulse profile up to the fourth harmonic was necessary. To do 
that we fitted each normalised pulse profile using the following 
expression: 



f(4>) = 1 + ^ fl « sin(2n7T(</> - (f>„)), 



(2) 



n=l 



where «i, «2, a 3 and 04 are the sinusoidal semi-amplitudes (here- 
after fractional amplitudes) of the first, second, third and fourth 
harmonics respectively and <f>i, $2, 4>3 and $4 the corresponding 
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Figure 2. In this figure a folded pulse profile is reported. For clarity two 
spin cycles are plotted. The continuous line is the best-fit using first, 
second and third harmonics. We also reported the single contribution 
to the profile of the first harmonic (dotted line), second harmonic (dot- 
dashed line), third harmonic (bi-dot dashed line) and fourth harmonic 
(dashed line). 



phases. We rejected the pulse phase delays for which the follow- 
ing two conditions were not both satisfied: i) the signal is not 
detected at least at a 3cr confidence level (Leahy et al. 1983); ii) 
the best fit fractional amplitude had to be at least at 3cr from zero 

(OLildCLi > 3). 

We first tried to fit the pulse phase delays with a polynomial 
to describe the pulse phase delays long term fluctuations plus 
the usual formula <p or b{t), describing the pulse phase residuals 
due to differential corrections to the initial orbital parameter es- 
timates(Deeter et al. 1981, see e.g. Riggio et al. 2007). 

We tried to describe the phase fluctuations using a poly- 
nomial up to the 9 th degree to fit obtaining a %r = 
6.1 1(1834.6/300), which is formally unacceptable. 

Due to the presence of these phase fluctuations we decided 
to apply the timing technique described in Riggio et al. (2007) 
in order to separate the orbital modulation from the phase fluc- 
tuations, to obtain a better estimates for the orbital parameters. 

Following Riggio et al. (2007) we define the pulse phase dif- 
ferences as: 



^<Porb(ti) = (porb(ti+l) ~ <Porb(ti), 



(3) 



We excluded all the points for which f, + i - f, > P or b in order to 
optimise the filter efficiency (Riggio et al. 2007). 

We iterated the process until convergence. In the last itera- 
tion (see Figure 3) the^ = 1.43(398/279), which is much more 
better than the previous approach. The orbital ephemeris best-fit 
results are reported in Table 1, where the errors have been mul- 
tiplied by the factor -Jx^ (Bevington & Robinson 2003). 

The values obtained are in perfect agreement with the ones 
reported by Papitto et al. (2009) from an analysis of a XMM- 
Newton observation and by Riggio et al. (2009), obtained with 
RXTE data but on a shorter time interval. The obtained uncer- 
tainties on the orbital parameters are a factor two smaller than 
the best orbital solution previously given by Riggio et al. (2009). 
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Figure 3. In the top panel of this figure the pulse phase delays (filled 
squares) and the corresponding pulse phase delay differences (filled cir- 
cles) on the data corrected with the orbital ephemeris obtained from the 
fit of the frequency Doppler shifts are reported. Each pulse phase point 
is obtained by fitting the folded pulse profile over about 1500 s of data. 
It is evident that there is no trace of such fluctuation in the pulse phase 
delay differences, giving a striking confirmation of the goodness of the 
method. In the bottom panel are reported the residuals (in <r units) with 
respect to the best fit orbital solution derived from the pulse phase delay 
differences using the timing technique described in the text. 



Note also that the orbital period we derive here has a relative un- 
certainty as small as 0.03%. 

In order to fit the pulse phase delays with a physically mean- 
ingful torque which takes into account the decreasing X-ray flux 
(and hence M, in the hypothesis that the X-ray flux is a good 
tracer of the mass accretion rate onto the NS), we need to de- 
scribe the flux evolution during the outburst, i.e. the bolometric 
flux at the peak of the outburst and the shape of the light curve. 

In order to estimate the bolometric flux of the source at 
the peak of the outburst, we considered the spectrum col- 
lected by the PCA in Standard 2 Mode (129 energy channels 
recorded at a time resolution of 16s) during Obs. 94041-01-01- 
04, which started on 55089.283 MJD and had a total exposure 
of 10480s. We considered only data taken by the PCU2 in or- 
der to avoid cross-calibration problem between PCUs. We se- 
lected only events detected in its top Xenon layer to maximise 
the signal to noise ratio(Jahoda et al. 2006). The background was 
modelled using the bright source model, that is appropriate for 
sources emitting > 40 c s _1 PCU -1 . We used the latest version 
( 1 1 .7) of the PCA response matrix generator, restricting to the 3- 
50 keV band and adding a systematic error of 0.5% to spectral 
counts 1 . Similarly to what was observed from other AMPs, the 
X-ray emission of IGR J1751 1-3057 is dominated by a power 
law-like emission extending to high (> 50 keV) energies. To 
model it we considered a simple thermal Comptonisation model 
(nthcomp, Zdziarski et al. 1996; Zycki et al. 1999), fixing the 



1 N. Shaposhnikov, K. Jahoda, C. B. Markwardt, 2009, 
http : / /www . universe . nasa . gov/xrays/programs/rxte/ 
pca/doc/rmf /pcarmf- 11.7/ 
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temperature of the hot electrons to 100 keV. Given the poor cov- 
erage of the PCA at low energies, we have also constrained the 
absorption column to nH = 1 x 10 22 cirT 2 , as it is suggested by 
an XMM-Newton observation performed during the same out- 
burst (Papitto et al. 2010). A 6.6 keV emission line was also 
added to model residuals in the iron range, though such a fea- 
ture is probably due to the contamination of the Galactic ridge 
in the field of view of the PCA (Markwardt et al. 2009b). The 
reduced chi square of the fit is good (y 2 = 64.3/72). The unab- 
sorbed flux we detect in the 3-50 keV band is 1.14(1) x 10~ 9 
erg cirT 2 s _1 , that extrapolated in the 0.5-200 keV band gives 
an estimate of ~ 2 x 10~ 9 erg cirT 2 s _1 . Under the hypothe- 
sis that emission is isotropic this corresponds to a bolometric 
luminosity of L x ~ 1.5 x 10 37 dg erg s _1 , where ds is the dis- 
tance to the source in units of 8 kpc. We note here that an upper 
limit of 10.6 kpc on the source distance was set by Bozzo et al. 
(2009), by imposing that the burst peak luminosity does not ex- 
ceed the Eddington limit, while Altamirano et al. (2010), from 
the analysis of the type-I bursts observed by RXTE and Swift, 
found an upper limit of 6.9 kpc. Papitto et al. (2010), on the ba- 
sis of the spectral analysis of XMM-Newton data, gave also a 
lower limit of ^ 7 kpc, although derived under some assump- 
tions. Moreover, as the source is only a few degrees away from 
the galactic centre and its X-ray emission is not heavily absorbed 
(nH ~ 10 22 cm 2 ), it is highly probable that the distance does 
not exceed 8 kpc. Assuming that L x = eGMM/R with e =* 1, 
we eventually deduce a peak mass accretion rate of the order 
of ^ 1.5 x 10" 9 dg M yr _1 , which is the estimate we use in 
the following to compare the dynamical estimates of M from the 
timing analysis. 

To describe the light curve shape, we chose to fit it with a 
piecewise linear function composed of three segments, as shown 
in Figure 1 . We choose to fit only the first 23 days of data since 
the subsequent data are affected by the concomitant XTE J 1 75 1- 
305 outburst. We modelled each of the three intervals with a 
function c,(0 = c,(l - (t - r,-)/T,-), T\ < t < T i+ i, where c* is 
the count rate at t = T, and t, the linear decay timescale for 
the i-th piece. For the sake of simplicity we wrote this piecewise 
function as cof(t), where cq is the count rate at the peak. The 
best-fit result is reported in Figure 1 . 



2.2. Timing analysis 

The spin frequency evolution in AMPs is thought to be driven 
by the accretion process. Matter falling from the accretion disk 
onto the NS transfers its angular momentum to the NS, which 
is spun-up to millisecond spin periods. But, as was evident from 
the first attempts (Ghosh et al. 1977), the magnetic field - ac- 
cretion disk interaction can exert a negative torque onto the NS, 
spinning it down. This is called the threaded disk model. Due 
to the complexity of the problem, the details of the NS mag- 
netosphere - disk interaction are still not well understood. In 
Literature three examples of AMPs which spin-down while ac- 
creting are reported (Galloway et al. 2002; Burderi et al. 2006; 
Papitto et al. 2008). 

As already observed in other two AMPs (see e.g. Burderi 
et al. 2007 for SAX J1808.4-3658 and Riggio et al. 2008 for 
XTE J 1807-294), also for this source the first harmonic is dom- 
inated by fluctuations and then unusable for our scope. It should 
be noted that the two AMPs cited above (SAX J 1808. 4-3658 and 
XTE J 1807-294) show a second harmonic with a more regular 
behaviour. 



An alternative interpretation of the pulse frequency deriva- 
tives was given by Hartman et al. (2008), who suggested that 
the red timing noise affecting the pulse phase delays can mimics 
a spin frequency derivative. Patruno et al. (2009) try to demon- 
strate that the pulse phase delays are correlated with the X-ray 
flux, rather than the genuine spin evolution of the source, due to 
motion of the hot-spot related to the flux. Unfortunately this cor- 
relation is not clear, even in the sign, in all AMPs and differs, in 
the same sources (see e.g. SAX J1808.4-3658 and XTE J1807- 
294), for each harmonic component, as noted by Patruno et al. 
(2009). We tested this hypothesis using the method as described 
in Patruno et al. (2009), adopting a constant spin frequency 
model to derive the pulse phase residuals. In the best-fit which 
maximise the linear correlation between phase residuals and flux 
we obtained^ 2 ^ 23 (4625 / 201 d.o.f.), indicating that, for this 
source, the pulse phase residuals of the first harmonic cannot be 
ascribed to flux variation and/or fluctuations. 

In the following, we will work under widely accepted hy- 
pothesis that the pulse frequency is the NS spin frequency. We 
will analyse the pulse phase delays and apply to them a disk 
threading model to derive the M and compare it with the value 
obtained from the spectral analysis of the same data. 

We epoch folded data on time intervals of about 3.0 ks (one 
pulse profile per data file) and 32 phase bins. A third harmonic 
is detected during all the outburst. A fourth harmonic was also 
detected in the first 15 days of the outburst. In Figure 4, 5, 6 and 
7 the first, second, third and fourth harmonic pulse phase delays 
are reported, respectively. The fractional amplitudes, as defined 
in Eq. 2, were corrected to take into account the instrumental 
background (~ 11 counts s _1 PCU 1 in the 2-25 keV energy 
band) and the background due to the presence of the galactic 
ridge in the field of view of RXTE (Markwardt et al. 2009b). To 
estimate this supplementary background we used the observa- 
tions of the AMP XTE J1751-305 when both sources went to 
quiescence, in particular the observations from MJD 551 15.400 
to 55126.745, with a total exposure of 33.7 ks. Due to the low 
count rate, we used the faint background model, obtaining a 
count-rate of ~ 7 counts s _1 PCU 1 in the same energy band. 
In Figure 8 the fractional amplitudes for the four harmonics are 
reported. 

The mean spin frequency reported in Table 1 is obtained by 
fitting the first harmonics pulse phase delays with a constant spin 
frequency model. The value obtained in this way for the spin fre- 
quency is v = 244.83395 15569(24) Hz. However, the systematic 
effects due to the uncertainty of 0.6" (Nowak et al. 2009) on the 
source position brings the mean spin frequency error to 7 x 10" 8 
Hz (see Burderi et al. 2007). 

2.2.1 . First harmonic 

We started fitting the first harmonic pulse phase delays with a 
constant spin frequency derivative in order to have an estimate of 
the mean spin frequency derivative. From the fit we obtained v = 
-7.0(9) x 10" 14 Hz s-' with a X 2 r = 18.51(3869/209), clearly 
unacceptable. 

In order to improve this result we used the threaded 
disk model applied to the AMPs by Rappaport et al. (2004). 
According to this model, the net torque acting on an AMP is: 

2 

r(f) = 2nIHt) = M(t) y/GMR c - ^ , (4) 

yK c 

where I is NS moment of inertia, M is the mass accretion rate, 
M is the mass of the NS, fi is the magnetic dipole moment of the 
NS and R c the co-rotation radius. 
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Figure 4. In the top panel we show the first harmonic pulse phase de- 
lays and the best-fitting curves, considering a constant spin frequency 
derivative (dashed line) and the disk threading model proposed by 
Rappaport et al. (2004). In the bottom panel we show the residuals (in 
cr units) of the first harmonic with respect to the best fitting model given 
by the threaded disk. 
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Figure 5. In the top panel of this figure the second harmonic pulse phase 
delays are shown. The dashed and the continuous lines are the second 
harmonic best-fit curves considering a constant spin frequency deriva- 
tive and the physical model considering the material torque proportional 
to the flux, respectively. In the bottom panel the residuals (in cr units) of 
the second harmonic with respect the constant spin frequency derivative 
are shown. 



To apply this model to our data we need an expression for 
M(t). We then made the hypothesis that the bolometric luminos- 
ity L(t) is a good tracer of M. In the hypothesis that the spectral 
variation during the outburst is not significant it is possible to 
consider L(t) proportional to the background subtracted count 
rate. We can then write M(t) = M m . ix f(t), where M max is the 
maximum accretion rate in correspondence of the flux peak and 
f(t) is the functional form of the count rate previously derived. 

An expression for v as function of M can easily be derived 
from Eq. 4: 



m 2/3pV3-Mr™ 

1.427 -^-^f(t)- 

'45 



■5.232 



/4 



mhsP\ 



xl(T 14 Hz s" 



where m is the NS mass in units of M Q , M™ 

-1-10 



(5) 

is the maximum 

mass accretion rate in units of 10~ lu M Q y _1 , P_3 the spin period 
in units of 1CT 3 s and I45 the NS moment of inertia in 10 45 gr 
cm 2 . In this work we adopted the FPS (Friedman Pandharipande 
Skyrme, see Friedman & Pandharipande 1981; Pandharipande 
& Ravenhall 1989) equation of state for which, fixing NS mass 
to M = 1.4M , we obtain a radius of R NS = 1.14 x 10 6 cm and 
a moment of inertia / = 1.29 x 10 45 gr cm 2 . 

The pulse phase delays formula used for the fit is obtained 
doubly integrating the expression 5 with respect to the time (see 
e.g. Burderi et al. 2007). 

From the fit of the first harmonic we obtained a magnetic 
dipole strength of fi = 1.64(7) x 10 27 G cm 3 and M max = 
5.7(6) x 10~ 9 M G y" 1 , with a X 2 = 12.36(2570/208), still unac- 
ceptable. The results are reported in Figure 4. Such a big x 2 ls 
clearly due to fluctuations in the pulse phase delays, as it is pos- 
sible to see in the best-fit residuals reported in Figure 4 (bottom 
panel). As foretold, we unsuccessfully try to interpret it with 
the model suggested by Patruno et al. (2009). The origin of such 
fluctuations still remains unexplained. 



2.2.2. Second harmonic 

Considering the second harmonic less affected by phase noise 
(Burderi et al. 2006; Riggio et al. 2008), we repeated the fit- 
ting procedure using these phase delays. In the constant spin fre- 
quency derivative case we obtained v = 1.45(16) x 10~ 13 Hz s _1 
with a. x 2 r — 1.74(238.5/137). The best-fit curve is reported in 
Figure 5. We adopted the same disk threading expression used 
for the first harmonic to describe the second harmonic. In this 
case we had to fix fi = since fi and M max strongly correlate. 
In this case, therefore, the derived value of M max has to be con- 
sidered as a lower limit to the mass accretion rate at the peak 
of the outburst, since value of the magnetic moment higher than 
zero will give a higher value for M max . The best-fit results are 
reported in Table 1. In this case the best-fit value of M max = 
0.92(10) x 10~ 9 M year 1 , with a X 2 = 1.70(232.8/137), is in 
good agreement with our estimate from the bolometric flux for a 
source distance of about 6.3(3) kpc. 

2.2.3. Third harmonic 

Since in this source the third harmonic is significantly de- 
tected in nearly the whole outburst we proceeded with the same 
method used for the first and second harmonic. From the fit 
with the constant spin frequency derivative we obtained v = 
244.83395151(6) Hz and v = 4.8(1.4) x 10~ 14 Hz s" 1 with a 
X 2 = 1.58(159.8/101). The best-fit curve is reported in Figure 
6. We applied the threaded disk model to describe the third har- 
monic and, as already done for the second harmonic, we fixed 
/j. — since /i and M max correlate in the fit. The obtained best-fit 
values are M max = 3.1(9) x 10~ 10 M year 1 for the peak mass 
accretion rate and v = 244.83395150(6) Hz with a reduced^ 2 
of^- 2 = 1.57(159.0/101). The third harmonic shows a spin-up, 
although not highly significant (> 3cr c.l.). The peak mass accre- 
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Figure 6. In the top panel of this figure the third harmonic pulse phase 
delays and its best-fit curves considering a constant spin frequency 
derivative (dashed line) and the physical model considering the material 
torque proportional to the flux (continuous line) are reported. In the bot- 
tom panel we show the residuals in units of a with respect the constant 
spin frequency derivative model. 
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Figure 7. In the top panel of this figure the fourth harmonic pulse phase 
delays and its best-fit curves considering a constant spin frequency 
derivative (dashed line) and the physical model considering the material 
torque proportional to the flux (continuous line) are reported. In the bot- 
tom panel we show the residuals in units of cr with respect the constant 
spin frequency derivative model. 



tion rate deduced is quite low in comparison to the one obtained 
from the bolometric flux for a source distance of 6.3(3) kpc. In 
fact the source distance is of 3.6(5) kpc. 

2.2.4. Fourth harmonic 

The fourth harmonic is sporadically detected in the first 15 days 
of the outburst, with a fractional amplitude of ~ 1%. From the 
fit with the constant spin frequency derivative we obtained |v| < 
2.7) x 1(T 13 Hz s- 1 (2cr c.l.) with a* 2 . = 4.03(44.28/11). The 
best-fit curve is reported in Figure 7. The obtained best-fit values 
are M max < 1.6 x 10~ 9 M year -1 (2cr c.l.) for the peak mass 
accretion rate with a reduced^ 2 of x 2 r = 4.09(45/1 1). 

It should be noted that, as already done in Burderi et al. 
(2007) and Riggio et al. (2008), we taken into account the effect 
of the source position uncertainty on the obtained values of v and 
v, adopting the same method described in Riggio et al. (2008). 
In particular, in the case of constant spin frequency derivative 
the uncertainties are of Av = 6.1 x 10~ 8 Hz on the frequency 
and Av = 0.72 x 10~ 14 Hz s _1 on the spin frequency derivative, 
while in the case of the physical model the uncertainties are of 
Av = 6.1 x 10~ 8 Hz on the frequency and AM = 4.3 x 10 -11 
Mo y _1 on the peak accretion rate. 

2.3. Spectral variability of the pulse profile 

As already done by Papitto et al. (2010) analysing the XMM- 
Newton observation of the same outburst, we analysed the en- 
ergy dependence of phase and amplitude of the three harmonic 
components. In order to have a good statistics we epoch folded 
data from 55089.233 MJD to 55092.5379 MJD, excluding the 
three data files where a type-I burst was present, with an expo- 
sure of 81 ks and a coverage of 28%. We choose this interval 




5 10 15 20 

Time (MJD - 55088.0) 

Figure 8. The fractional amplitudes of the four harmonics are reported. 
Points in correspondence of the type-I bursts are not plotted due to dif- 
ficult estimate of the persistent flux. 



because the pulse phase delays have a quite stable linear trend 
with an average first harmonic phase scattering of 0.005. We 
consider only data from PCU2. We evaluated also the contri- 
bution to the background given by the galactic bulge, using the 
data described in the previous section, in the working hypothesis 
that the changes due to the different observation data and slightly 
different instrument pointing are small (see Papitto et al. 2010, 
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Figure 9. In the left panel the phase lags for the first three harmonics are reported. The phase lags are measured with respect to the maxima of 
each component, as done in Papitto et al. (2010). In the right panel the fractional amplitude for the three harmonics are reported. 



for a detailed discussion). For these reasons small systematics 
could be present and the confidence intervals could be underes- 
timated. The fractional amplitude of the first harmonic (Fig. 9, 
right panel) increases from 17.5(3)% at 2.5 keV to 26.7(3)% at 
6 keV. It remains roughly constant around 23% up to 12 keV, 
where the fractional amplitude is 25.4(3)%. Above 12 keV there 
is a steady decline of the fractional amplitude, although up to the 
energy band 25-60 keV is still well detectable with a fractional 
amplitude of 5.2(3)%. The second and third harmonic fractional 
amplitudes show a behaviour similar to the first harmonic one, 
increasing with energy with a maximum (around 10-20 keV) of 
3.7(4)% and 2.2(4)%, respectively. We also observe a decline 
of the fractional amplitude above 20 keV for these harmonics. 
Fractional amplitudes as a function of photon energies were not 
clearly detected for the fourth harmonic, probably due to the 
long time of integration and the smearing caused by the observed 
fourth harmonic phase fluctuations of ~ 0.05 (see Fig. 7), which 
is a considerable fraction of the fourth harmonic period. 

The first harmonic show phase lags (reported in Fig. 9, left 
panel). As described in Papitto et al. (2010), there is a steady 
decrease of the phase lag up to ~ 10 keV, where it is clearly 
visible a break. Beyond 10 keV the pulse phase lag still decrease, 
but with a smaller rate. 

The second harmonic shows a different trend with respect to 
the first harmonic. It reaches the maximum lag around 10 keV, it 
shows no time lag at 25 keV, but then we observe a sudden jump 
of ~ -450/is from 25 to 50 keV. 

The third harmonic phase lags are roughly constant in all the 
energy band. It should be noted that in the energy band 25.7- 
59.8 keV the phase lags for the second and the third coincide. 

2.4. type-l bursts timing 

We performed a timing analysis of all the 10 type-I burst present 
in the observation, with the goal of studying the pulse profile 



evolution during the burst. In Figure 10 we show the results for 
the second and fifth type-I burst, the best sampled among the 
bursts present in the RXTE observation. The starting date of these 
bursts is 55089.721 MJD (TDB) and 55094.619 MJD (TDB), 
and the decay time is 8.0(1) seconds and 8.5(1) seconds, respec- 
tively. We divided each burst in chunks holding (roughly) the 
same number of events so that in each folded profile a fractional 
amplitude of ~ 20% can be easily detectable. We folded each 
chunk using 8 phase bins and performed an harmonic decompo- 
sition using only the first harmonic. The results of this analysis 
are reported in Figure 10. During both bursts, the pulse phase 
delays remain stable and, with the exception of the very first 
seconds, locked to the pulse phase delays during the persistent 
emission (see Figure 10, mid panel). The fractional amplitude 
behaviour is even more interesting because remains, within the 
errors, quite constant during both bursts and locked to the per- 
sistent value (Figure 10, bottom panel). A detailed spectral and 
temporal analysis of all these type-I bursts was recently reported 
by Altamirano et al. (2010), although with different techniques. 



3. Discussion 

Among the AMPs known, IGR J1751 1-3057 shows several pe- 
culiarities. The source shows the highest persistent first har- 
monic fractional amplitude with a peak value (background cor- 
rected) at the beginning of the outburst of about 23% which 
linearly decrease to 17% at the end of the outburst (see Figure 
8), while the highest ever observed fractional amplitude was ob- 
served by Patruno et al. (2010) in XTE J 1807-294. 

The pulse shape is complex, showing, on integration times 
of 3 ks, a second harmonic, a third harmonic and sporadically 
a fourth harmonic with fractional amplitudes of 2.5%, 1.6% and 
1%, respectively (see Figure 8). A third harmonic as strong as the 
second harmonic could be, following Poutanen & Beloborodov 
(2006), an indication that both hot spots are visible with the sec- 
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ondary spot only partially visible, since for a single spot the third 
harmonic should be much smaller than the second harmonic 
(02/03 > 5, Poutanen & Beloborodov 2006), suggesting inter- 
mediate values for the inclination angle. In the AMPs, Hartman 
et al. (2008) reports of a sporadically detectable third harmonic 
in SAX J1808.4-3658, while Patruno et al. (2010) detected spo- 
radically a third and a fourth harmonic in XTE J 1807-294. 

The fact that in IGR J1751 1-3057 the third harmonic is vis- 
ible for nearly all the outburst with a total of 105 detections over 
216 folded pulse profiles makes this source peculiar. 

Moreover, while the fractional amplitude of the first har- 
monic clearly shows a steady decrease with the flux, such a 
decrease is less evident in the second and third harmonics, for 
which the fractional amplitude remains more stable when the X- 
ray flux decreases (see Figure 8). 

However, the most interesting and puzzling result is the 
different behaviour of the phase delays of the four harmonic 
components. In particular, similarly to other AMPs, such as 
SAX J1808.4-3658 and XTE J1807-294 (see e.g. Burderi et al. 
2007; Hartman et al. 2008, 2009; Riggio et al. 2008; Patruno 
et al. 2010), the first harmonic shows clear phase fluctuations 
(see Figure 4). A discussion over the results based on the anal- 
ysis of the first harmonic is strongly affected by these phase 
fluctuations. We can only do some hypotheses on the nature of 
such fluctuations. The amplitude of these fluctuations, ~ 0.05 in 
phase units, could corresponds to hot spot movements of ~ 18 
degrees, already seen in numerical simulations (see Romanova 
et al. 2004; Bachetti et al. 2010), although on time-scales of 
fraction of seconds, while the observed time-scales are of few 
days. We also tried to interpret these phase fluctuations with the 
model suggested by Patruno et al. (2010). We adopted a con- 
stant spin frequency model and a linear relation between flux and 
pulse phase residuals. The best correlation parameter set gave a 
X 2 r - 23, clearly unacceptable. A model to describe and correctly 
interpret these fluctuations is still missing. 

However, under the working hypothesis that some exchange 
of angular momentum between the NS and the accreting matter 
has to occur during X-ray outbursts, when the accretion rate is 
at its maximum, we try to interpret the behaviour of the phase 
delays of the second and third harmonic, which appear to be less 
affected by phase fluctuations. The fact that the phase delays de- 
rived from the second harmonic appear more stable than those 
derived from the first harmonic has been observed, for instance, 
during the 2002 outburst from SAX J1808.4-3658. Burderi et al. 
(2007) showed that, while the first harmonic phase delays clearly 
show a phase shift at days 14 from the beginning of the 2002 out- 
burst, a similar phase shift was not present in the phase delays 
derived from the second harmonic. A similar behaviour was also 
observed by Riggio et al. (2008) for the AMP XTE J 1807-294, 
which went onto outburst just once in 2003 in the RXTE era. In 
both these cases, the interpretation of the phase delays derived 
from the second harmonic in terms of accretion torques provided 
reasonable spin frequency derivatives (and inferred mass accre- 
tion rates onto the NS), although a different interpretation was 
given for both sources by Hartman et al. (2009); Patruno et al. 
(2009). Naturally, these results have to be taken with great cau- 
tion, since it is not excluded that phase fluctuations can still af- 
fect phase delays derived from the second (or higher) harmonic. 

In the case of IGR J 1751 1-3057, again, the second harmonic 
shows a more regular behaviour with respect to the first har- 
monic and suggests a spin-up of the NS (see Figure 5). From the 
fit with our simplified torque model we obtain an M estimate, fix- 
ing n = 0, of M max = 0.92(10) x 10~ 9 M year 1 , which would 
be compatible with the observed X-ray flux from the source if 



we put the source at a distance of 6.3 kpc. This could indirectly 
suggest that the second harmonic component is a better tracer of 
the spin frequency evolution, even if a physical model that ex- 
plains all the phenomenology observed in all the AMPs class is 
still missing. Results obtained on the second harmonic give, as 
in the case of the AMPs SAX J1808.4-3658 (Burderi et al. 2006; 
Hartman et al. 2009) and XTE J1807-294 (Riggio et al. 2008), 
reasonable values for the physical parameters of the system. 

However, if we look at the third harmonic we find that it 
shows a spin-up although it is not highly significant (> 3<x) and 
smaller than the value inferred from the second harmonic. The 
lower limit to the accretion rate M max = 3.1(9) x 10" 10 M year -1 
gives a distance of 3.6(5) kpc, at ~ 5<x from the value inferred 
from the second harmonic. This value is not compatible with the 
second harmonic one. The lack of literature and observations of 
the third harmonic behaviour in other AMPs does not allow to 
make a comparison with other cases, leaving the question open. 

A fourth harmonic was sporadically detected during the out- 
burst. We tried to fit the pulse phase delays with the same models 
adopted for the other harmonics but the result is inconclusive. It 
should be noted that is affected by a phase fluctuations (~ 0.02) 
comparable with the first harmonic's fluctuations. It is probable 
that integrating on long time-scales, the large fluctuations make 
its detection unfeasible. 

We performed a high resolution timing analysis of all the 
type-I bursts present in this RXTE observation. The results are 
very similar for all the bursts, and we show the burst 2 and 5 
(see Figure 10) since these are the best sampled in the RXTE 
observation. The first harmonic phase delays appear to rise in 
correspondence of the fast rising phase of the burst, suggesting a 
frequency drift during the first few seconds since the burst onset 
(Altamirano et al. 2010), and to return to the phase value during 
the persistent emission during the burst decay, giving some evi- 
dence that the burst probably starts not far from the hot spot in 
the polar cap. The fractional amplitude in each burst remains in- 
stead locked (within the errors) to the persistent emission value 
during all the burst, suggesting the quite surprising fact that the 
temperature gradient does not vary during the burst. An interpre- 
tation of a such behaviour is beyond the scope of this work. 

From the orbital ephemeris reported in Table 1 the pulsar 
mass function is f x = 1 .070854(21) x 10~ 3 M . From this value 
of mass function we can derive a minimum mass for the compan- 
ion star of 0. 14 M , considering an inclination angle of 90° and a 
NS mass of 1.4 M . With such a minimum mass, the companion 
star of IGR J1751 1-3057 is one of the more massive compan- 
ion stars among the AMPs, together with XTE J 1807-294 and 
SAX J1748.9-2021 (Altamirano et al. 2008). Using the relation 
R RLl = 1.2 x lO 10 ™ 1 ^^ 3 cm(Paczynski 1971), where m 2 , .i 
is the companion mass in O.1M units and Pih is the orbital pe- 
riod in two hours units, we obtain for the companion's Roche 
lobe radius a value of 0.248 R . Such value is larger than what 
is expected for a low mass main sequence star (see Chabrier & 
Baraffe 2000, 5 Gyr track), for which the corresponding radius 
is about 0.15 R . It can be shown that the contact condition be- 
tween Roche lobe and companion star pose a firm lower limit to 
the inclination of the system of ~ 20 degrees, corresponding to 
a companion mass of ~ O.45M . For smaller inclination angles 
the companion star would overfill its Roche lobe. This obviously 
excludes that the companion star could be a white dwarf or an 
helium-core star, while strongly suggests that the companion star 
is a main sequence star, possibly bloated as a consequence of 
its evolutionary history (Podsiadlowski et al. 2002). The nature 
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Figure 10. Two type-I burst timing analysis is reported. For each burst in the top panel the PCU2 count rate of the type-I burst present in this 
observation, in the mid panel the pulse phase delays and in the bottom panel the fractional amplitude are reported. 



of the companion star is thoroughly discussed in Papitto et al. 
(2010). 
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